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Abstract 

The conductance in Coulomb blockade quantum dots exhibits sharp peaks 
whose spacings fluctuate with the number of electrons. We derive the 
temperature-dependence of these fluctuations in the statistical regime and 
compare with recent experimental results. The scrambling due to Coulomb 
interactions of the single-particle spectrum with the addition of an electron 
to the dot is shown to affect the temperature-dependence of the peak spacing 
fluctuations. Spectral scrambling also leads to saturation in the temperature 
dependence of the peak-to-peak correlator, in agreement with recent experi- 
mental results. The signatures of scrambling are derived using discrete Gaus- 
sian processes, which generalize the Gaussian ensembles of random matrices 
to systems that depend on a discrete parameter - in this case, the number of 
electrons in the dot. 
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A quantum dot is a sub-micron-sized conducting device containing up to several thousand 
electrons. In closed dots, the coupling between the dot and the leads is weak and the charge 
on the dot is quantized [Tj. The addition of an electron into the dot requires a charging 
energy of Eq = e 2 /C (where C is the capacitance of the dot). This charging energy can 
be compensated by varying the gate voltage V g , leading to Coulomb blockade oscillations of 
the conductance versus V g . In the quantum regime (i.e., for temperatures below the mean 
level spacing A), conductance occurs by resonant tunneling, and sharp conductance peaks 
are observed function of V g . 

Dots can be fabricated with little disorder such that the electron dynamics in the dot 
is ballistic. Larger dots are often characterized by irregular shape, resulting in chaotic 
classical dynamics of the electrons. Such dots are expected to exhibit universal mesoscopic 
fluctuations which are the signature of quantum chaos. In particular, the distributions 
of the conductance peak heights in Coulomb blockade quantum dots at T <C A)) were 
predicted to be universal, depending only on the underlying space-time symmetries || . The 
measured distributions were found to agree well with theory. The statistics of the 
peak heights at finite temperatures (T ~ A) were also derived recently using random matrix 
theory (RMT) |5|]. The measured distributions || become narrower and less asymmetric with 
increasing temperature, in qualitative agreement with theory, although significant deviations 
were observed at higher temperatures, presumably due to dephasing. 

Another quantity whose statistics was recently studied both experimentally and theoreti- 
cally is the spacing between successive conductance peaks. In the simple constant interaction 
(CI) model (where Coulomb interactions are included only as an average charging energy) 
and for T A, a (shifted) Wigner-Dyson peak spacing distribution is expected, but the 
observed distributions are Gaussian-like f?HiO|. This has been explained as an interaction 
effect by numerical diagonalization of a small Anderson model with Coulomb interactions 
|?|,|1T||. The temperature dependence of the peak spacing statistics was also measured re- 
cently ||. In this paper we use the finite-temperature theory plus RMT to study the peak 
spacing fluctuations at temperatures T ~ A. We find a rapid decrease of the fluctuations 
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above T/A ~ 0.5, in agreement with the experimental results. 

Interaction effects beyond the charging energy were not included in the finite temperature 
theory. They can be treated in a single-particle framework within a mean-field approximation 
(e.g., Hartree-Fock). Due to charge rearrangement, we expect the spectrum to change or 
"scramble" upon the addition of an electron into the dot |T2|,|5|. The effects of scrambling 
on the statistics can be described by a random matrix model that depends on a discrete 
parameter: the number of electrons on the dot. The theory of discrete Gaussian processes 
|j~3| can then be used to analyze the finite temperature statistics of peak spacings and peak 



heights for various degrees of scrambling. A rescaled parametric distance controls how fast 
the spectrum is changing, and will be referred to as the scrambling parameter. It was 
shown that spectral scrambling can lead to nearly Gaussian peak spacing distributions at 



low temperatures [Q. In this paper we derive two main signatures of a changing spectrum 
in the finite temperature statistics: the less rapid decrease of the spacing fluctuations with 
temperature for T/A > 0.5, and the saturation of the number of correlated peaks at higher 
temperatures |15|| . The first effect has not been experimentally observed while the second 
has been qualitatively suggested and observed in Ref. ||. We also derive a simple expression 
for the scrambling parameter in terms of the dot's properties. 

At T ~ A, several resonances contribute to each conductance peak [16| 



G(T, E F ) = j = w x (T, E F )g x , (1) 

A 

where g is the dimensionless conductance expressed as a thermal average over individual level 
conductances g\ = 2T~ l T l > T r x / (T l x + T r x ). The thermal weight w\(T,Ep) of a level A (for 
T < E c ) is given by w x = Af(AF M -E F )(n x )Jl-f (E x - E F )]. AF M = F{M)-F{M-1) 
where Fjj is the canonical free energy of M non-interacting electrons on the dot, Ep = 
E F + eaVg — (Af — 1/2)E C is an effective Fermi energy (a is the ratio between the plunger 
gate to dot capacitance and the total capacitance), and (n\) is the canonical occupation 
of a single-particle level A. Both the canonical free energy and occupation are calculated 
exactly using particle- number projection ||. In the statistical theory, the eigenvalues E\ 



and wavefunctions ip\ fluctuate according to the corresponding Gaussian random matrix 
ensemble. The fluctuations of the partial widths T\ are calculated by relating the widths to 
the eigenfunctions across the dot-lead interfaces. 

Eq. (|1|) was used in Ref. to calculate the conductance peak distributions by full 
RMT simulations, and an approximate analytic expression was derived in the limit where 
spectral fluctuations are ignored. The finite temperature formulation can also be used to 
calculate the temperature-dependence of the peak spacing distributions. Unlike the peak 
height statistics, the peak spacing statistics are sensitive to fluctuations of both the spectrum 
and the wavefunctions, and full RMT simulations are required. The location of the A/"-th 
peak is determined by finding the value of Ep for which the conductance ([3]) is maximal. 
Statistics of peak spacings are collected from different successive peaks as well as from 
different realizations of the dot's Hamiltonian. The peak spacings exhibit less fluctuations at 
higher temperatures as is demonstrated in the top panel of Fig. [I], where the spacings A2 for 
a typical peak series calculated in one random matrix realization are shown at T/A = 0.5 and 
T/A = 2. While we do not expect to reproduce the observed functional form (i.e. Gaussian- 
like) of the peak spacing distribution using a fixed spectrum, it is still meaningful to study 
the temperature-dependence of the standard deviation of the spacings cr(A 2 ) (where A 2 = 
(A2 — (A2))/A). The results for the GUE statistics are shown in the bottom panel of Fig.[l] 
(solid line). The width shows a slight increase until about T/A ~ 0.5 and then decreases 
rapidly with increasing T/A. We compare the calculations with recent experimental results 
H in the presence of a magnetic field (circles). The calculations somewhat underestimate the 
experimental width but describe well the overall observed temperature dependence. Shown 
in the inset is the calculated ratio <7goe(A2)/ctgue(A2) which increases as a function of 
T/A. The experimental ratio of ~ 1.2 — 1.3 measured at T ~ 100 mK is consistent with our 
theoretical results. 

So far we have taken into account only an average value for the Coulomb interaction 
{M 2 Ec/2). Interaction effects can be described within a single-particle framework in the 
Hartree-Fock (HF) approximation [|l7j]. The peak spacing at T <C A can be expressed as a 



second order difference of the ground state energy of the dot as a function of particle number 
0. According to Koopmans' theorem, the change in the HF ground state energy when an 
electron is added is given by the HF single-particle energy of the added electron 

MM+l) _ p{N~) _ P (JV+1) / 9 x 
°HF °HF ~ ^Af+l i y^J 

where S^ F is the HF ground-state energy of the dot with i electrons, and Ej is the energy 
of j-th single-particle state for a dot with i electrons. The theorem is valid in the limit 
when the single-particle eigenf unctions are independent of the number of electrons, and is 
expected to hold for large dots and for interactions that are not too strong. Its validity 



in numerical models of quantum dots was recently studied JT^j. The spacing A 2 (Af + 1) 
between the Af-th and Af + 1 peak is then given by 

A 2 (Af + 1) = - = AE^ + AE N . (3) 

AE^ +1 ^ = EjH[^ — E^^ is the level spacing for a fixed number of electrons (Af+ 1), and 
AEj^f = E^ +1 ^ — E^ is the energy variation of the AA-th level when the Af + 1 electron is 
added to the dot. 

The HF Hamiltonian of the dot depends on the number of electrons, and in the following 
we denote by H(xj^) the Hamiltonian for Af electrons (in this notation E^ = E\(xjf)). 
For a dot whose single-electron dynamics is chaotic we shall assume that H(xj^) is a discrete 
Gaussian process, i.e., a discrete sequence of correlated Gaussian ensembles of a given sym- 
metry class ||13|| . Such a process can be embedded in a continuous process H{x) where the 



Hamiltonian depends on a continuous parameter x [18, IS]. A parametric dependence that 



originates in the dot's deformation as a function of V g was used to explain the Gaussian-like 



shape of the peak spacing distribution flTJj]. Here the parametric dependence is assumed 



to be mainly due to interaction effects, as recent experimental results indicate |6|;[TI|. We 
further assume that Ax = x/v+i ~ X N is approximately independent of Af. According to 
the theory of Gaussian processes, the parametric statistics are universal upon scaling of 



the parameter by the rms level velocity ||20|| . We denote by Ax the variation of the scaled 



parameter between successive number of electrons. The rms level velocity depends on the 
symmetry class, and in the following the parameter is always scaled by the GUE rms level 
velocity. 



A simple Gaussian process is given by [fig, IS 



H{x) = cosrcf/i + smxH 2 , (4) 

where Hi and H2 are uncorrelated Gaussian random matrices. For each realization of 
the Gaussian process we calculate the single-particle spectrum E\( function of the 

parameter. At T ~ A, Eq. (§) does not hold. Instead we determine the spacing from 
the location of successive peaks: the M-th peak is determined using the levels E\{x^f) and 
wavefunctions ip\{xj^) as explained before, while the M + 1 peak is determined similarly 
but using a different spectrum E\{x^f + \) and eigenf unctions i^x( x M+i)- Consequently, the 
spacing A2 depends now on both T/A and Ax. 

Fig. |2| shows the standard deviation of the GUE peak spacing distribution <t(A 2 )/A as 
a function of T/A (on a log-log scale) for several values of the "scrambling" parameter Ax. 
As for the Ax = case, we observe a decrease above T/A ~ 0.5, except that the decrease is 
more moderate. It would be interesting to see whether this signature of spectral scrambling 
can be observed experimentally. 

Another quantity that is sensitive to spectral scrambling is the peak-to-peak correlator, 
which is characterized by its full width at half maximum (FWHM), i.e., the number of 
correlated peaks n c . For a constant single-particle spectrum, n c is found to increase linearly 
with temperature since the number of levels contributing to a given peak is ~ T/A. However, 
for a changing spectrum the number of correlated peaks is expected to saturate (as a function 
of temperature) at a value ~ m that measures the number of electrons needed to scramble 
the spectrum completely. This effect was observed in the experiment |J. It can be calculated 
quantitatively by using the Gaussian process (f|) |I5]| . For each value of Ax and T/A, the 
peak-to-peak correlations are determined universally. The top left panel of Fig. ^ shows the 
calculated peak-to-peak correlator c(n) at several temperatures for a constant spectrum (i.e., 



Ax = 0), where the correlator width is seen to increase with temperature. The right inset 
shows the correlator c(n) for the same temperatures but for a changing spectrum Ax = 0.5, 
where its width is seen to saturate at higher temperatures. The bottom panel of Fig. [3] shows 
the number of correlated peaks n c versus T/A for several values of Ax. n c saturates at a 
smaller value of m as Ax gets larger, i.e., when the single-particle spectrum scrambles faster. 
To further illustrate this effect, we show in Fig. f| the peak height fluctuation g — g for a 
series of peaks using one realization of the GP @ at T/A = 0.5 and 2. For a fixed spectrum 
(Ax = 0) we observe a significant increase of the correlations at the higher temperature (left 
panels), while for a changing spectrum (Ax = 0.5) the correlations do not change much with 
temperature (right panels). 

Experimentally it was found that, in a smaller dot, n c saturates at a smaller value m 
||. This suggests faster scrambling (i.e., larger Ax) in the smaller dot. We can derive an 
expression for the scrambling parameter Ax in terms of the dot's properties by relating 



the parametric approach to the microscopic HF-RPA approach of |I2| . In the latter the 
fluctuations in the variation AEj^- of the A/"-th level upon the addition of an electron (see 
Eq. (|3|)) are dominated by charge that is pushed to the surface, and are estimated to be 
a 2 (AEtf) ~ ft~ 1 A 2 /g, where g is the dimensionless Thouless conductance. In the parametric 
approach we find, using perturbation theory fl9 |, <j 2 (AEj^) = 2ft~ 1 A 2 (Ax) 2 . Comparing 



the two expression for a 2 (AE^) we obtain 

Ax oc g- 1 ' 2 ~ M~ 1/A . (5) 

The symmetry class parameter ft drops out in this relation. Indeed we expect Ax to depend 
only on the dot's properties irrespective of the presence or absence of a magnetic field. 
Relation (^) is valid in the regime where a 2 (AEj^-) is linear in (Ax) 2 , i.e., for Ax < 0.3 (see 
Fig. 2 of Ref. [0]). The second estimate in (|5|) is for a ballistic dot where g ~ Af 1 ^ 2 . Relation 
(H) confirms that Ax is larger for the dot with a smaller number of electrons. Notice the 
qualitative similarity between the theoretical results of Fig. |3] and the experimental results 
in Fig. 2 of Ref. §. 



In conclusion, we used the statistical theory of Coulomb blockade quantum dots to cal- 
culate the temperature-dependence of the peak spacing fluctuations. Statistical scrambling 
of the spectrum upon adding an electron to the dot affects the temperature-dependence of 
both the peak spacing fluctuations and the peak-to-peak correlator. 

This work was supported in part by the Department of Energy grant No. DE-FG-0291- 
ER-40608. 
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FIG. 1. top panel: A typical peak spacing A2 = (A2 — (A2))/A sequence obtained from one 
RMT realization of the dot's Hamiltonian at T/A = 0.5 and T/A = 2. Notice the reduced 
fluctuations at the higher temperature. Bottom panel: The standard deviation of the peak spacing 
cr(A2) for the GUE statistics as a function of T/A using a log-log scale. The solid line is the 
theoretical result and the circles are the experimental results of || for two dot configurations. The 
inset is the calculated ratio o"goe(A2)/o"gue(A2) as a function of T/A. 
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FIG. 3. Spectral scrambling and the finite temperature peak-to-peak correlation. Top panels: 
the peak-to-peak correlator c(n) as a function of distance n between peaks for several temperatures 
and for a fixed spectrum Ax = (top left panel) and a changing spectrum with Ax = 0.5 (top 
right panel). Bottom panel: the number of correlated peaks n c (defined as the full width at half 
maximum of the correlator c(n)) as a function of T/A for several values of the scrambling parameter 
Ax = 0,0.1,0.2,0.5 and 1. Notice the saturation of n c at smaller values for larger values of Ax 
(i.e., smaller dots; see Eq. (|5|). Compare this figure to the experimental figure 2 of Ref. Q. 
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FIG. 4. Spectral scrambling and peak height fluctuations. Shown is the peak height fluctuation 
g — g for a series of peaks generated from one realization of a Gaussian process. Left: T/A = 0.5 
and T/A = 2 for a fixed spectrum (Ax = 0). Right: T/A = 0.5 and T/A = 2 for a changing 
spectrum (Ax = 0.5) 
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